Dynamic crossover scaling in polymer solutions 
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The crossover region in the phase diagram of polymer solutions, in the regime above the overlap 
concentration, is explored by Brownian Dynamics simulations, to map out the universal crossover 
scaling functions for the gyration radius and the single-chain diffusion constant. Scaling considera- 
tions, our simulation results, and recently reported data on the polymer contribution to the viscosity 
obtained from rheological measurements on DNA systems, support the assumption that there are 
simple relations between these functions, such that they can be inferred from one another. 
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The basic physics of systems of linear flexible un- 
charged macromolecules in solution, in a wide range of 
concentrations and solvent qualities, can be considered as 
reasonably well understood, at least as far as the stan- 
dard scaling laws are concerned, which govern the static 
and near-equilibrium dynamic properties of the solution 
[ll-Q- The fairly rich structure of the generic phase dia- 
gram of this system [1,13, shown in Fig. [TJ is a result of 
the interplay between molecule-molecule overlap on the 
one hand, and solvent quality effects on the other. These 
factors determine whether excluded-volume (EV) inter- 
actions between monomers are important or not. For 
highly dilute systems, where overlap is unimportant, we 
may characterize the effective strength of the EV inter- 
action (or quality of the solvent) by the dimensionless 
parameter z* = (e/^T)(l — 0/T), where ks is Boltz- 
mann's constant, T is the temperature, 6 is the Theta 
temperature at which an isolated infinitely long chain 
would collapse, and e is the repulsive energy (of enthalpic 
origin) with which two monomers interact. For small z*, 
EV effects are unimportant, and the chain conformations 
obey random- walk (RW) statistics (R ~ WV 1 / 2 , where R 
is the size of the chain, N its number of monomers and b 
the monomer size) , while for large z* the chain is swollen 
and obeys the statistics of a self-avoiding walk (SAW), 
R ~ bN v where v « 0.6 is the Flory exponent. Upon 
increasing the concentration c (number of monomers per 
unit volume) beyond the overlap value c*, Flory screen- 
ing sets in, such that for such systems again R ex TV 1 / 2 , 
even in the case of a good solvent. This interplay is 
well described by the "blob" concept p], 0, 0] , where the 
thermal blob size £t characterizes the length scale (in 
real space) up to which EV is unimportant due to an 
insufficiently large value of z*, while the overlap (or con- 
centration) blob size £ c is the typical length scale above 
which chain overlap induces Flory screening. A semidi- 
lute solution (regime C in Fig. [T]) is characterized by 
b <C £t <C £ c <^ R, i- e. by only a finite window of length 
scales between £x and £ c where SAW statistics apply. 
Upon further increasing c, the system enters the concen- 
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FIG. 1. (Color online) Phase diagram of a polymer solution 
in the plane monomer concentration c and excluded-volume 
strength z* , as predicted by the standard blob picture. The 
(-) regime occurs for small z* <C iV -1//2 , while z* ^> iV _1//2 in 
the good solvent regime. The overlap concentration is given 
by c*. For concentrations above c**, excluded-volume inter- 
actions are completely screened. For derivations of relations 
governing the crossover boundaries, see Supplemental Mate- 
rial [ti- 



trated regime c ^> c**, where this window has shrunk 
to zero. The lines drawn in Fig. [I] thus do not indicate 
sharp phase transitions but rather smooth crossovers. 

These crossovers in the static conformational proper- 
ties are accompanied by a dynamic crossover [lHj]- For 
c <C c*, hydrodynamic interactions (HI) are important, 
such that the dynamics is characterized by the scaling 
laws of the Zimm model (most notably D ~ kBT/(r] s R), 
where D is the single-chain diffusion constant, and ?7 S 
the solvent viscosity). For c ^> c*, HI are screened, es- 
sentially as a result of chain-chain collisions which tend 
to randomize the momentum propagation in the system. 
The fact that HI in semidilute solutions are unscreened 
for short times [6] (up to the Zimm relaxation time of an 
overlap blob, which is the typical waiting time for chain- 
chain collisions to occur) is of no relevance for the present 
paper, since we are interested in the long-time dynamics 
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FIG. 2. (Color online) Same as Fig. [T] but with c/c* and 
z = z^N 1 ^ 2 as variables. In this representation the loci of 
the crossovers do not depend on the chain length N. The 
points where the present simulations and recently reported 
experiments [12] have been performed are indicated as well. 



only. Therefore here the scaling laws of the Rouse model 
(e. g. D oc N~ r ) apply. We consider only systems where 
N and c are still small enough that entanglements play no 
essential role, such that a further crossover to reptation 
dynamics [l|-[4| does not need to be taken into account. 
HI screening is governed by the same length scale £ c as 
EV screening jajTzl- 

Applying the blob picture to the various regimes of 
Fig. [TJ one finds straightforwardly the scaling results for 
£t> £c> R, D, and the polymer contribution to the viscos- 
ity 7?p, listed in Table I of the Supplemental Material [Ej], 
and we refer the reader to the literature [4] for further 
details. 

In the past, these scaling laws (which are asymptotic 
results applicable only sufficiently far away from the 
crossover regions), as well as the concentration-driven 
crossovers, have been subjected to a wealth of carefully 
executed experiments [4j, [8| and computer simulations 
0, 0]. However, to the best of our knowledge, there 
have been only few investigations that have tried to 
study the double crossover driven by both solvent quality 
and concentration, in particular when dynamic proper- 
ties are concerned [lOj. The present Letter attempts to 
fill that gap. We are particularly interested in the uni- 
versal crossover scaling functions , which describe the 
behavior in the asymptotic limit N — >• oo. 

The first step in analyzing universal crossover behavior 
is the introduction of suitably scaled dimensionless vari- 
ables (or ratios). We therefore re-parameterize the scal- 
ing laws by replacing the EV strength z* with the "sol- 
vent quality" parameter z = z^N 1 / 2 (note that this is the 
standard crossover scaling variable used to describe the 
O transition), and the monomer concentration c with the 



ratio c/c*. The results of this transformation are shown 
in Fig. [2] and listed in Table HI For R, D and r] p we in- 
troduce dimensionless ratios as indicated in the table. In 
this representation the dependence on the chain length 
TV has been completely absorbed in the N dependence of 
c* and z, and the only remaining variables are c/c* and 
z. It is therefore natural to generalize these results to the 
relations R/R* = 4>r(c/ c* , z\ D/D* = <Pd{c/ e* , z), and 
Vp/Vp = ^r?( c / c ^5 z \ where i?*, D*, rj* denote the values 
of R, D, r] p at c = c*, and (/>£>, and ^ are universal 
crossover scaling functions defined on the whole plane of 
Fig. [21 and which, up to numerical prefactors, are known 
deep in the asymptotic regimes (these are just the re- 
sults listed in Table I]) , but whose behavior needs to be 
calculated or measured near the crossover lines. 

An interesting observation can be made when look- 
ing at the last two lines of Table H In all regimes above 
the overlap concentration, we find that the functions 
R/R*, D/D*, and rj p /rj* are not independent, but can 
be calculated as soon as one of them is known. It 
is then natural to assume that the same property also 
holds in the crossover regimes, i.e. that the relations 
(c/c*) 1 / 4 ^^^ 4 = const, and (c/c*) 1 / 3 ^^ 1 ^ 6 = const, 
hold in the whole plane c/c* ^> 1 — or, in other words, 
that there is essentially only a single crossover scaling 
function in the problem. 

In what follows, we present results for <j)R and c^d ob- 
tained by careful computer simulations. It turns out 
that the convergence to the universal behavior is quite 
rapid, and can be treated fairly reliably by extrapola- 
tions [13[. In addition, we use recently reported data 
on (j)^ (obtained from rheological measurements on DNA 
systems) [12[, to support the hypothesis of the existence 
of a single crossover scaling function. Due to the univer- 
sality of our results, they should be broadly applicable to 
a large class of experimental systems. 

A linear bead-spring chain model is used to repre- 
sent polymer molecules, with each polymer chain coarse- 
grained into a sequence of N beads connected by massless 
springs. An ensemble of N c such bead-spring chains, en- 
closed in a cubic cell of edge length L and immersed in an 
incompressible Newtonian solvent, is used to model the 
polymer solution. The monomer concentration of beads 
per cell is c = N to t/V, where N tot = N x N c is the total 
number of beads, and V = L 3 is the volume of the sim- 
ulation cell. A Brownian dynamics algorithm [3] which 
accounts for HI between the beads, and which can simu- 
late both 6 and good solvent conditions, is used here to 
obtain the stochastic time evolution of bead positions. In 
particular, the pair- wise EV interactions between beads 
v and \i in a chain are represented by a narrow Gaussian 
potential TH, 15[ 
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TABLE I. Various normalized quantities of the polymer system (thermal blob size £t, overlap blob size £ c , polymer radius 
R, single-chain diffusion constant D, polymer part of the viscosity ?7 P ) in the regimes indicated in Fig. [21 in terms of c/c* and 
z. Blob sizes are not indicated where they are irrelevant. Numerical prefactors of order unity have been ignored. R* , D* , 77*. 
denote the values of R, D, n p at c = c*. For derivations, see Supplemental Material |5|. 



where d* is a non-dimensional parameter that mea- 
sures the range of excluded volume interaction, b is 
the monomer size, and r UfI is the vector connecting 
monomers v and \i. Simulations to explore the solvent 
quality crossover have been carried out by keeping the 
value of z constant at four different values: {0, 0.7, 1.7, 3}, 
with z = corresponding to 6-solutions. At each 
fixed value of z, values of c/c* ranging from 0.1 to 4 
have been used to sample both the dilute and semidilute 
regimes. The overlap concentration c* is calculated from 
c* = N/ [(47r/3)i?p ] > where R g o is the radius of gyration 
of a chain in the dilute limit, which is computed apriori 
by running single chain BD simulations. The box size L 
is selected to ensure L > 2R e (where R e is the end-to- 
end distance of a chain), in order to prevent a chain from 
wrapping over itself. For this purpose, the value of R e 
at any value of c/c* is estimated from the blob scaling 
law R e = R e o ( c / c *) - ( 2l/-1 )/( 6i/-2 ) where R e o is the end- 
to-end distance of a chain computed in the dilute limit. 
Typical simulations consist of an equilibration run for 
approximately one relaxation time followed by a produc- 
tion run of 10 to 60 relaxation times; here, the relaxation 
time r was estimated via r = R^/(6D), where R g is the 
concentration-dependent radius of gyration. Moreover, 
data were obtained by averaging over 30 — 70 indepen- 
dent runs. 

Asymptotic results at any point (c/c*,z), independent 
of simulation parameters, are obtained by accumulating 
finite chain data from BD simulations carried out at fixed 
values of these variables, and subsequently extrapolating 
to the limit N — >• 00. At finite values of TV, the results still 
depend on the parameters d* [see Eq. (pQ)], and on the hy- 
drodynamic interaction parameter h* = a/ ^/ttQ^bT / H), 
where a is the bead radius and H is the spring constant. 
The parameter <i* is truly irrelevant in the long-chain 
limit, i.e. asymptotic conformation statistics are only 
governed by c/c* and z [13|- Conversely, h* is a variable 
similar to z*, where the relevant scaling variable govern- 
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FIG. 3. (Color online) Universal predictions of the ratio 
D/D z at the two state points (z, c/c*) = (0.7, 3) and (1.7, 2), 
where D is the single-chain diffusion constant at the state 
point under consideration, while Dz is the c = value at the 
same z. Data accumulated for finite chain lengths N (sym- 
bols) for two different values of the hydrodynamic interaction 
parameter (h* = 0.15 and 0.28) extrapolate to a unique value 
in the limit N — >• oo. The limited range of N values leads to 
a conservative estimate of the error of the extrapolated value. 



ing the crossover from free-draining to non-draining be- 
havior is the draining parameter h = h*\/N [13j |. In the 
present study we focus on the non-draining limit h — » oo, 
which is obtained by TV — >• oo at constant h*. This limit is 
independent of h*, as demonstrated in Fig. [3] for D/Dz, 
which is the ratio of the single-chain diffusion constant 
at a finite concentration to its value in the dilute limit. 
Data for chains of lengths TV = 6 to N = 20 have been 
extrapolated to TV — > oo. 

The resulting crossover scaling functions (/>r(c/c*,z) 
and cVd(c/c*, z) are shown in Figs. [H where the mean 
polymer size R is R g . The data are limited to the 
crossover region where in principle none of the power laws 
of TableUhold. Since the observation regime is small, it is 
nevertheless possible to fit a power law to the data (with- 
out deeper theoretical justification). This gives rise to a 
convenient description in terms of an effective empirical 
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FIG. 4. (Color online) Concentration dependence of the 
crossover scaling functions <j)R and <j)D in the semidilute 
regime, for different values of solvent quality z, obtained 
by Brownian dynamics simulations. The effective exponents 
zy e ff have been determined by fitting power laws to the data 
cj)R oc (c/c*)~ p and requiring p = (2z/ e ff — l)/(6^ e ff — 2), accord- 
ing to Table U in regime C. Using these values, lines are drawn 
according to (\>r> oc (c/c*)~ q with q = (1 — ^ e ff)/(3^ e ff — 1)- <\>r> 
data for z — are omitted because of large statistical errors. 




FIG. 5. (Color online) Demonstration that the functions 
R/R* , D/D* , and ?7p/?7p are not independent, but can be cal- 
culated as soon as one of them is known. In the upper panel, 
z — data are omitted because of large statistical errors. 



exponent (for values see Figs. |4]), which is obtained 
by adjusting the value of v in the scaling laws of regime 
C of Table U to the observed slope. Since the crossover 
scaling functions are related to each other (see below), 
each crossover yields the same value for ^ e ff- Moreover, 
consistency with the asymptotic laws requires that ^ e ff 
varies between 0.5 and 0.6. A similar crossover scaling 
behavior is obtained for ^(c/c* , z) from measurements 



of rj p for semidilute DNA solutions across a range of N 
c, and z* 

The validity of the expectation from scaling consid- 
erations that the combination (c/c*) 1 / 4 (j)R (j) 1 ^ = 1 in 
the entire semidilute regime, while it is proportional to 
(c/c*) 1//4 in the dilute limit, and similarly that the com- 



1 in the entire semidilute 



bination (c/c*) 1 / 3 </)r (j>n 1//6 
regime and (c/c*) 1 / 6 in the dilute limit, is established 
in Figs. [5] In the case of the latter combination, while 
values of </)r have been obtained from the present simula- 
tions, values of ^ v are from the data reported in Ref. Il2. 
Clearly, the knowledge of a single two-argument scal- 
ing function is sufficient to describe the entire double 
crossover regime in polymer solutions. 
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BLOB SCALING IN THE CONCENTRATION 
AND SOLVENT QUALITY CROSSOVER 
REGIMES 

In regimes B and C of the phase diagram (see Fig. 1 in 
the main paper), where solvent quality effects play a role 
in determining the static and near-equilibrium dynamic 
properties of the solution, the length scale correspond- 
ing to thermal blobs is relevant. The size of the thermal 
blob £t can be estimated as follows. Two monomers are 
assumed to interact with an energy e (1 — 0/T), T > 0, 
e > — i. e. the interaction is repulsive. On small 
scales, the interaction is not strong enough to perturb 
the random-walk (RW) nature of the chain. The chain 
notices that there is repulsion only when the total repul- 
sive energy is of order /c#T. This defines the thermal 
blob size, which we assume corresponds to n monomers. 
Since the number of contacts in a RW is of order n 1 / 2 in 
three spatial dimensions, the condition to find n is 



,1/2 



k B T. 



(1) 



In the present work, the pair-wise excluded- volume (EV) 
interactions between monomers v and \i in a chain are 
represented by a narrow Gaussian potential, 
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where, d* is a non-dimensional parameter that mea- 
sures the range of excluded volume interaction, b is 
the monomer size, and r vyb is the vector connecting 
monomers v and \i. The parameter z* is the non- 
dimensional strength of the pair-wise excluded-volume 
interactions. It follows that, 
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k B T 

Equation (pQ) then implies, 

71= (z*)~ 2 . 

The corresponding blob size is consequently 

£ T = bn 1/2 = b(z*)- 1 . 
Using the definition of the solvent quality parameter, 



(3) 



(4) 



(5) 



where, N is the number of monomers in a chain, the 
non-dimensional blob size in terms of scaled variables is, 



bN 1 / 2 



(7) 



In regimes B and C, on scales larger than £t excluded 
volume is important. In regime B, where the monomer 
concentration c < c*, the entire chain conformation is 
that of a self-avoiding walk (SAW) on length scales larger 
than However, in regime C, where c > c*, Flory 
screening sets in at yet larger length scales leading to the 
existence of concentration blobs. Thermal blobs within 
the concentration blob obey SAW statistics. If £ c is the 
size of such a concentration blob, then its size is given by 



(8) 



where, v is the Flory exponent, while m is the number of 
thermal blobs within the concentration blob. The size £ c 
is found by the knowledge that the solution is homoge- 
neous on scales larger than £ c , and that the concentration 
blobs are space filling. The total number of monomers 
within £ c is nm, therefore 



nm = cQ. 



(9) 



Using Eqs. (jl]), J5J), and J5]), it is straight forward to show 
that 

m=(z*c- 1 b- 3 ) 1/i3u ~ 1) . (10) 
As a result, using Eqs. © and (JTUJ) in Eq. (jSJ 

2u-l v 



£ c = 6(2*) 3«/-l(c& d ) 3*-l. 



(11) 



In order to represent the size of the concentration blob 
in terms of scaled variables, it is necessary to calculate 
the overlap concentration c*, which separates regimes B 
and C. The overlap concentration can be found from the 
condition that at c = c*, the concentration blob contains 
just all monomers of the chain, N. From Eqs. (j9]), (j8j) 
and (j5j), since nm = N at c = c* 



N = c*b 6 z 



(12) 



Substituting for m from Eq. ([10]), and solving for c* , leads 
to 



z = z+N 1 ' 2 , 



(6) 



c* = &- 3 iV-( 3i/ - 1 ^*- 3(2i/ - 1 ). 



(13) 
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When solved for z*, this leads to the equation of the line 
that separates regimes B and C in the unsealed variables 
phase diagram (Fig. 1 in the main paper). 

It is convenient to represent cb 3 as (c/c*)(c*6 3 ) when 
converting from unsealed to scaled variables. One can 
then show from Eqs. ([TT|) . (|6j) and ([T3]) that the non- 
dimensional size of the concentration blob in terms of 
scaled variables is 



bN 1 / 2 



(?) 



(14) 



Upon systematically increasing the concentration c, 
the concentration blob shrinks until at a threshold con- 
centration c = c** , £ c = £r, and there is then no longer 
any length scale where SAW statistics apply. The locus 
of points c** separates regime C from regimes D and E 
in the unsealed variables phase diagram. It is straight 
forward to calculate c** by equating Eqs. (|5j) and JTTJ 

c** = &" 3 **. (15) 

In terms of scaled variables, it is similarly straight for- 
ward to show by equating Eqs. ([7]) and ([T4]) that 



„2(3i/-l) 



(16) 



In regimes D and E, where only RW statistics apply, 
the concentration blob size is given by 



- &„V2 



Zc = bg 



(17) 



where, g is the number of monomers in a concentration 
blob. Since the solution is homogeneous on all length 
scales larger than £ c , and the concentration blobs are 
space filling, 
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(6^1/2)3 



which implies 



and as a result 



9 = (cb 3 )- 2 



Hc = b(cb 3 )- 



(18) 



(19) 



(20) 



The expression for c* given by Eq. ([T3]) is only valid away 
from the regime, i.e., for z* > TV -1 / 2 . On the other 
hand, in the O regime, since at c = c* the concentra- 
tion blob contains all the monomers in the chain, which 
implies £ c = 67V 1 / 2 , Eq. ([20]) leads to 



= b- 3 N-^ 2 . 



(21) 



This difference in the expressions for c* makes it neces- 
sary to distinguish between regimes D (where EV effects 
are important) and E (where EV effects are negligible) 
when representing the size of the concentration blob (and 



indeed all other observables) in terms of scaled variables. 
In regime D, Eqs. ([20]) . (fT3l) and (|6]) imply 



-(?)"'■ 



,3(2i/-l) 



bN l/2 Vc x 

while in regime E, Eqs. ([20]) and ([2T]) imply 



bN 1 / 2 



\c* 



(22) 



(23) 



Once the blob scaling laws for £x are known in regimes 
B and C, and those for £ c are known in regimes C to E, 
the scaling laws for all other observables in these regimes 
can be derived. Even though these laws have been de- 
rived and discussed in detail previously (H-Q, here, as ex- 
amples, we derive the scaling laws for the static chain size 
R, and the single-chain diffusion coefficient D, in order 
to represent them in terms of the notation used in this 
work. The scaling laws for regime A are not discussed 
since the conventional notation in this regime is followed 
here. Once the expressions for R and D are known, the 
scaling laws for r] p , the polymer contribution to viscosity, 
can be derived from the relation r] p ~ ksT (c/N) r, where 
r is the longest relaxation time of the macromolecule, 
r - R 2 /D. 



Regime B 

Since the chain obeys SAW statistics on length scales 
larger than 



R = ir ( - 

n 



Using Eqs. (JU and (j5]), leads to 

R = bN v {z*) 2v - 1 . 



(24) 



(25) 



The scaling relation in terms of scaled variables can be 
obtained by making use of Eq. (jlj) 



R 



bN 1 / 2 



(26) 



Zimm dynamics are relevant on all length scales. 
From the Nernst-Einstein relation, D = kBT/rj s R, and 
Eq. (|23 



D 



r] s bN u (z*) 



(27) 



where, ?7 S is the solvent viscosity. Defining a non- 
dimensional diffusion coefficient D / {ksT /riJjN 1 ^ 2 )^ it 
follows from Eqs. ([6]) that it is given in terms of scaled 
variables by 



D 



(k B T/ Vs bN^) 



-(2I/-1) 



(28) 
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TABLE I. Various quantities of the polymer system (thermal blob size £t, overlap blob size £ c , polymer radius R, single-chain 
diffusion constant D, polymer part of the viscosity rj p ) in the regimes indicated in Fig. 1 of the main paper, as a function 
of monomer size 6, chain length N, excluded- volume interaction strength z*, thermal energy fc^T, solvent viscosity rj s , and 
monomer concentration c, within the framework of scaling theory. Blob sizes are not indicated in cases where they are irrelevant. 
Numerical prefactors of order unity have been ignored. The scaling laws are valid in the asymptotic regimes sufficiently far 
away from the crossover boundaries. 



Regime C 

Since RW statistics are obeyed on length scales larger 
than £ c 

1/2 



R = tc 



N 



nm 



(29) 



where, N/nm is the number of concentration blobs in a 
chain. Using Eqs. (@]), JlOj and ([TT]) . leads to 



1 2v-l 



R = bN 1 ' 2 [^(cfr 3 )- 1 ] 23 -- 1 . 



(30) 



The scaling relation in terms of scaled variables is ob- 
tained by using Eqs. (|6]) and ([13} 
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bNV< 
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1 2v-l 
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(31) 



Rouse dynamics are obeyed on length scales larger 
than £ c . Therefore 

k B T 



D rj s £ c (N/nm) ' 
Using Eqs. (g]), (TO) and dTTJ) , leads to 



(32) 



D 



k B T , . 



VsbN 



(z*y^v-i( c tf)-3v-i. (33) 



The non-dimensional diffusivity in terms of scaled vari- 
ables follows from using Eqs. (j6j) and (fT3|) in Eq. ([33]) 
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ik B T/r] s bN^) 
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l-v 



(34) 



Regime D 

In regime D, RW statistics are obeyed on all length 
scales, both within the concentration blob, and by the 
concentration blobs themselves. As a result 



(35) 



where Eqs. ([T9]) and ([20]) have been used for simpli- 
fication. The non-dimensional chain size is trivially 
R/ibN 1 / 2 ) = 1. 

Rouse dynamics are obeyed on length scales larger 
than £ c . Using Eqs. ([32]) (with g = nm), and Eqs. ([T9j) 
and (EH leads to 



D 



r!sbN {C0) ' 



(36) 



The non-dimensional diffusivity in terms of scaled vari- 
ables follows from Eqs. ([6]) and ([T3]) 



D 



(37) 



{UbT /rjsbN 1 / 2 ) \c 



Regime E 

In regime E, the scaling relations for R and D in terms 
of unsealed variables are identical to those in regime D, 
since in both regimes, RW statistics are obeyed on all 
length scales, and Rouse dynamics apply on length scales 
larger than £ c . In terms of scaled variables, while the 
distinction between the two regimes is irrelevant for the 
static chain size, the difference in the expression for c* 
due to the presence and absence of EV effects, respec- 
tively, manifests itself as a difference in the expression for 
the non-dimensional diffusivity. Using Eqs. ([6]) and ([2T]h 
Eq. ([36]) reduces in terms of scaled variables to 



D 



ikBT/risbN 1 /*) 



\c* 



(38) 



Table U summarises all the unsealed equations derived 
in this section, including those for the polymer contribu- 
tion to viscosity. Table I in the main paper summarises 
the equations in terms of scaled variables. 
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ESTIMATION OF THE ERROR IN THE 
DIFFUSION COEFFICIENT D 

Error at finite N 

The long time self-diffusion coefficient is calculated 
from the mean-square displacement (MSD) of the cen- 
ter of mass of each chain [4]. Each stochastic trajectory 
in the simulation leads to a time series for the MSD. 
An ensemble average over all the trajectories then gives 
the time series of the mean of the MSD (denoted here 
as MSD avg ). Since the initial MSD avg data represents 
transient short time diffusivity, while the data at large 
times have large error bars (they are based on a smaller 
number of blocks in the block averaging method), we typ- 
ically discard the first 15 - 20% and last 20 - 30% of the 
data. This leads to a window of times At in the MSD avg 
data that can be fitted with a straight line. At a given 
value of TV, the size of Ar depends on the state point 
(z,c/c*). For a fixed set of parameters, the slopes of the 
lines fitted to each of the MSD trajectories in the ensem- 
ble, over the range of times Ar, is used here to obtain 
an ensemble of predicted diffusion coefficients. The mean 
diffusion coefficient D and the standard error of mean is 
then determined from this ensemble. 

Error in the extrapolated value at N —> oo 

At each state point (z,c/c*) the mean value of D and 
the error in the mean, for a set of finite size chains with 
N = 6,8,..., 20, is obtained as described above. The 
mean and error-of-mean of the ratio D/Dz is then ob- 
tained via single chain BD simulations to determine Dz 
for the same values of N and z. As described in the 
main paper, the asymptotic value of D/Dz in the long 
chain limit is obtained by plotting the data at finite N 
as a function of iV -1 / 2 , and extrapolating a straight line 
fit to the data, to the limit N — >> oo. The fitting of 
the straight line, and the estimation of the error in the 
extrapolated value is carried out here with the help of 
"Least-Squares Fitting" numerical routines provided in 
the GNU Scientific Library (GSL). 

Essentially, the GSL routine finds the least-squares fit 
by minimizing x 2 , the weighted sum of squared residuals, 
for the straight line model D/Dz = Co + c\N~ 1 ^ 2 . The 
weights are the inverse of the error at each data point. 
The fitting routine gsl_fit_wlinear returns the best-fit pa- 
rameters Co and c\ , along with a 2 x 2 covariance matrix 
that measures the statistical errors on Co and c\ result- 
ing from the errors in the data. The standard deviations 
of the best-fit parameters are then given by the square 
root of the corresponding diagonal elements of the co- 
variance matrix. Particularly conveniently, the routine 
gsl_fit_linear_est uses the best-fit coefficients cq, c\ and 
their estimated covariance to compute the fitted func- 



tion and its standard deviation at any desired point. By 
using gsl_fit_linear_est to find the value of the fitted func- 
tion and its error at N~ x l 2 = 0, we determine both the 
infinite chain length limit value of D/Dz and its associ- 
ated error. 

SCALING OF COMPUTATIONAL COST WITH 
CHAIN SIZE 

As illustrated in Fig. 3 of the main paper, asymptotic 
predictions in the long chain limit have been obtained 
for each state point (z, c/c*) by extrapolating finite chain 
data accumulated for chain lengths ranging from N = 6 
to N = 20. The use of this rather limited range of N 
values is necessitated by the computational cost of the 
Brownian dynamics simulation algorithm used here. As 
shown below, an estimate of the computational cost of 
the algorithm can be derived by using some simple scaling 
arguments, and the resulting expression can be verified 
by comparison with the CPU time required in the current 
simulations. 

The computational cost of carrying out an Euler inte- 
gration of the governing stochastic differential equation 
for a single time step has been shown in Ref. |4| to scale 
with system size as N£ ot , where x = 2.1. Recall that 
^V~tot = N x N c is the total number of beads in a cu- 
bic cell of edge length L, with N c being the number of 
bead-spring chains. Since typical simulations consist of 
runs extending over several relaxation times r, followed 
by averaging over many independent runs, it is necessary 
to find the dependence of r on N in order to find the 
scaling of the total CPU cost. This is done as follows for 
simulations carried out in the semidilute regime C. 

From Eq. (|3T|) . the requirement that L > 2R in order 
to prevent chains from wrapping over themselves leads 
to 

1 2v-l 

L~b(±>) ""-V- 1 ^/ 2 . (39) 

The concentration of monomers in the simulation box is 
c = (NN C )/L 3 . As a result, from Eq. ([39]) 

-(f)^) 11 ^^- 1 ^" 172 - (40) 

Using Eqs. ([T3]) and (j6]), the overlap concentration c* can 
be written in terms of the solvent quality z as 

c* ~ b- 3 z-^ 2v -^N- 1 / 2 . (41) 

It follows from Eqs. ([40]) and (|4T|) that N c is independent 
of z, and related to the scaled concentration through the 
relation 

Nc-^Y^. (42) 
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FIG. 1. (Color online) CPU time as a function of chain length 
for Brownian dynamics simulations carried out at the two 
state points (z, c/c*) = (0.7,3) and (1.7,2) for two different 
values of the hydrodynamic interaction parameter (h* = 0.15 
and 0.28), on a 156 SGI Altix XE 320 cluster. Symbols rep- 
resent the various values of chain length N used in the simu- 
lations, while the dashed line is drawn with a slope predicted 
by the scaling relation Eq. (|44[) . 

The number of chains in a simulation box is consequently 
constant when (c/c*) is maintained constant. 

The relaxation time of a macromolecule r ~ R 2 /D. 
From Eqs. (|3T|) and this implies 



2-3v 



TQ ( c / c *) 3^T z Wv-i) N 3/2 m 



(43) 



where To = rj s b 3 /ksT is the monomer relaxation time. 

The total computational cost of a single stochastic tra- 
jectory consequently scales as r (N C N) X . For fixed values 
of z and (c/c*) therefore 



and (1.7, 2), for two different values of the hydrodynamic 
interaction parameter /i*, on a 156 SGI Altix XE 320 
cluster. The simulations are identical to those used to 
display the ratio D/Dz in Fig. 3 of the main paper. It 
is immediately apparent that the prediction in Eq. ([44]) 
for the scaling of the total CPU time with chain size is 
obeyed closely by the present simulations, independent 
of solvent quality and scaled concentration. 

It is clearly desirable to add data for longer chains in 
order to improve the accuracy of the asymptotic value ob- 
tained by extrapolation to the long chain limit. Unfortu- 
nately, the additional CPU cost this would entail (which 
can be estimated from the scaling of computational cost 
with chain size shown in Fig. [I]) makes this highly infeasi- 
ble. For instance, simulating a chain with N = 30 would 
require roughly (9(10 4 ) CPU hours, while N = 100 would 
require roughly O(10 6 ) CPU hours! Indeed, by including 
simulations for a chain with N = 30, the CPU time for 
obtaining the static and dynamic properties of a semidi- 
lute solution at a single state point (z,c/c*) — (3,4) in 
the phase diagram (Fig. 2 of the main paper) is estimated 
to increase from its current value of approximately 7 x 10 4 
hours to roughly 3 x 10 5 hours. 
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hosted by University of Melbourne. 



Total CPU time per run 



3 i 



TV" 



-7.2 



(44) 



where in the last expression on the right-hand-side we 
have substituted the value x = 2.1 for the current algo- 
rithm, and used an exponent for N that enables a repre- 
sentation of the CPU cost as displayed in Fig. [U 

The various symbols in Fig. [T] represent simulations 
carried out at the two state points (z,c/c*) = (0.7,3) 
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